Can Cosmic Ray Catalysed Vacuum Decay 
Dominate Over Tunnelling? 

§^ ; Kari Enqvist^ and John McDonald^ 

5h . Department of Physics, P.O. Box 9, FIN-00014 University of Helsinki, Finland 



< 



February 1, 2008 



^ ■ Abstract 

I We consider the question of whether cosmic ray catalysed false vacuum decay 

^ ■ 

Q_i' can be phenomenologically more important than spontaneous decay via quan- 

Qh! tum tunnelling. We extend the zero bubble wall width Landau- WKB analysis 

' of catalysed false vacuum decay to include the leading order effects of finite 

- ^ , wall width and derive an expression for the thin- wall bubble action. Using 

' this we calculate the exponential suppression factor for the catalysed decay rate 

at the critical bubble energy, corresponding to the largest probability of catal- 
ysed decay. We show that, in general, cosmic ray catalysed decay is likely to 
be more important than spontaneous decay for sufficiently thin-walled bubbles 
(wall thickness less than about 30% of the initial bubble radius), but that spon- 
taneous decay will dominate for the case of thick-walled bubbles. Since any 
perturbative model with a cosmologically significant false vacuum decay rate 
will almost certainly produce thick-walled bubbles, we can conclude that cos- 
mic ray catalysed false vacuum decay will never dominate over tunnelling in 
imposing phenomenological constraints on perturbative particle physics models. 
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1 Introduction 



False vacuum decay ^ is known to play an important role in constraining the param- 
eter space of many models of particle physics. Examples include the Standard Model 
itself 0, which can have a metastable vacuum state for a range of the parameters 
which determine the Higgs potential, and the Minimal Supersymmetric extension of 
the Standard Model (MSSM) 0, which can have charge and colour breaking min- 
ima for a range of its parameters. In particular, the question of the stability of the 
electroweak vacuum with respect to decay to colour breaking minima is of particu- 
lar significance for the possibility of successful electroweak baryogenesis within the 
framework of the MSSM §. 

Usually it is assumed that the most important false vacuum decay mode is sponta- 
neous decay via quantum tunnelling to the true vacuum state |T| . However, it has been 
recognized that false vacuum decay could also be catalysed, in principle, by collisions 
of cosmic rays which produce a virtual scalar field with sufficient energy to nucleate 
a super-critical bubble of false vacuum p. Early attempts to estimate the rate of 
such catalysed events indicated that such processes might be more important than 
spontaneous decay for a range of the parameter space 0. However, these discussions 
were based upon somewhat uncertain physical reasoning and assumptions, making it 
difficult to assess their validity. Later, the problem was approached both numerically 
and analytically. The most recent numerical analysis, based on a lattice calculation of 
the decay rate, indicates that the exponential factor in catalysed decay is about 0.8 
times that in spontaneous decay [0. This calculation was done in the context of a 
— theory and so corresponds to the extreme case of a thick-walled bubble. 

In the present paper we will be concerned with the analytical treatment of catalysed 
false vacuum decay proposed by Voloshin based on the Landau- WKB method of 
calculating transitions between strongly different states in quantum mechanics 0. In 
this approach one writes a Lagrangian for a bubble of true vacuum as a function of 
its radius R, L{R, R), and then quantizes this action and computes the amphtude for 
the creation of a bubble of radius R from the initial false vacuum state by using the 
WKB method as developed by Landau. Our main goals in this paper will be to extend 
the Voloshin calculation, which was based on bubbles of vanishing wall thickness, to 



the case where the effects of bubble wall width are included and then to apply our 
results to the question of whether cosmic-ray catalysed false vacuum decay can be 
more important that spontaneous decay in imposing phenomenological constraints on 
particle physics models. Compared with the earlier attempts to address this question, 
we believe that our approach provides a much more firm and quantitative answer. 

The paper is organized as follows. In section 2 we review the application of the 
Landau- WKB method to false vacuum bubble formation in the presence of a virtual 
scalar field for the case of bubbles of zero wall width. In section 3 we derive an action 
for thin-walled bubbles including the leading order corrections due to finite wall width. 
In section 4 we apply the Landau- WKB method to this action and calculate the effect 
of finite wall width on the bubble formation probability. In section 5 we discuss 
cosmic ray catalysed and spontaneous false vacuum decay and the phenomenological 
implications of our results. In section 6 we present and discuss our conclusions. In 
the Appendix we review the original Landau discussion of the application of the WKB 
approximation to the calculation of matrix elements of strongly different states in 
quantum mechanics. 

2 Virtual Scalar Field Catalysis of False Vacuum 
Decay in the Zero- Width Limit 

In this section we set up our notation and review the application of the Landau- WKB 
method to the catalysis of false vacuum decay in the limit of vanishing bubble wall 
thickness [Q. We will consider a real scalar field with a potential U{(f)) of the form 

U{^) = ^i4>' - v'f - M^ct> ^ UM - M'4> , (1) 

which has two minima with an energy difference e = 2vM^. For further use we define 
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ft = ^V2Xv'^ . (2) 
o 3 

It can be shown that the classical Lagrangian for a zero-width bubble wall reads 



L{R, R) = -fl\/ l-R^ + iR^ , (3) 



where R = R{t) is the radius of the bubble. In the next section we derive and 
generahze this resuh for non-zero bubble width. The classical Lagrangian Eq. @ has 
one degree of freedom, which may be canonically quantized as for the case of ordinary 
quantum mechanics. In particular, the WKB method of Landau for the calculation 
of the transition matrix element between two "strongly different" states may be 
applied to Eq. (^. (We review the quantum mechanical Landau- WKB method in the 
Appendix). 

Landau showed that, up to a relatively unimportant prefactor, the transition matrix 
element of an operator / between two states ipi and 1^2 {E2 > Ei) is given by 



/12 ~ exp 



Re I r ^2m{U - Ei)dx - f ^2m{U - E2)d. 



(4) 



where oi and 02 are the classical turning points of the motion of the particles in the 
presence of a potential barrier U{x) and x* is the "transition point", corresponding to 
a point in the complex plane at which the derivative of the exponent with respect to 
X* is zero. This expression can be directly generalized to the matrix element of any 
operator between strongly different initial and final states H, 



(F(E2)|/|X(Ei)) ~exp [-lL]=exp 



Crq* rqy ^ 

i / p{q; Ei)dq + i i p{q; E2)dq 
Jqx Jq* I 



(5) 

where p are the canonical momenta, qx and qy are the classical turning points and g* 
is the transition point. In the following we will refer to li^ as the Landau integral. 

For the case of the bubble Lagrangian of Eq. (|^), the canonical momentum p/j is 
given by 

dL liR^R 
PR = —r= I > (6) 



and the Hamiltonian H = prR — L is given by 



H = -jtE= - ~eR' . (7) 



In order to evaluate the Landau integral, we must solve the Hamiltonian Eq. in 
order to find R as a function of R and E and so obtain as a function of R and E. 
This step turns out to be particularly easy for the zero-width case. We obtain 

R^ = J- (8) 



and 

Pn = {E'-f,'Rr/\ (9) 

where E = E + eR^. 

The amphtude < B{E)\0 >mrtuai for the transition from the false vacuum state 
to a bubble of energy E in the presence of a virtual scalar field is related to the 
matrix element of the scalar field < B{E)\(j){x)\0 > integrated over x. This can be 
seen by noting that a virtual scalar field is created by including a source term J(j) 
in the Lagrangian. Thus, representing the amplitude by a path integral, we see that 
switching on the source J will produce a change in the amplitude according to 

s_sim2: . jL| B(Em^)\o > . (10) 

In the effective theory of thin wall bubbles of radius R the non-trivial part of the field 
operator integrated over all space becomes 



J d^x{(j){x) + v) 



-R-". 11 



Thus the calculation of the matrix element of the scalar field operator reduces to the 
calculation of the matrix element of R^ in the effective theory. Thus in terms of the 
dimensionless variables x, e and ^, defined by 

R = x^, E = e^, e=4' (12) 

the exponential suppression of the amplitude for the creation of a super-critical bubble 
of energy E from the vacuum state (corresponding to a bubble of zero energy and 
radius) in the presence of a virtual scalar field is given by ||^ 



{B{E)\R^\0) ~exp (-^Re 



rx* fxiE) I 

j y x'^ — x^dx + J y x"^ — {x^ + eydx 



(13) 



The turning point x{E), corresponding to the radius of a bubble of energy E, is 
obtained by setting i? = in the Hamiltonian Eq. (^, which gives 

x{Ef-x{Ef = e. (14) 



We note that, whereas in the original Landau derivation the positive square root should 
be taken in the canonical momentum (see the Appendix), for the case of super-critical 



bubble formation from a sub-critical bubble initial state the directions of the incoming 
momenta of the initial and final bubble states are opposite, implying that the negative 
square root should be taken in the second integral. This is also necessary in order to 
obtain a non-trivial complex transition point. The transition point is then given by 
the solution of 



Vx^ - x6 = ^x^ - (x3 + e)2 . (15) 

Thus in the zero width case the transition point is given by x*{e) = {-e/2y/\ In- 
tegrating up to the root in the upper half-plane, Voloshin showed that the smallest 
exponential suppression factor, corresponding to the largest probability of catalysed 
vacuum decay, is obtained when the energy corresponds to the maximum energy for 
which a sub-critical bubble can exist i.e. the energy of the critical "sphaleron" bubble 
configuration at the top of the energy barrier seperating the sub-critical and super- 
critical bubble regions [Q. From H{R = 0) we see that this corresponds to a bubble 
of radius Xc = 2/3 and energy Cc = 4/27 ~ 0.15. The exponential suppression factor 
is then given by 

|(fi(E,))|i?3|0)r~exp(-F(Ec)) , (16) 

where 

^ 0.16 (17) 
F(0) ^ ' 

and F{0) corresponds to the bounce action for vacuum tunnelling |jl[]. For energies 
greater than the critical energy the additional energy will be dissipated via perturbative 
particle production (which is favoured over formation of a higher energy bubble), such 
that no significant increase in the vacuum decay probability is obtained for E > Ec 
i- 

In the above calculation the exponential suppresion factor is, in fact, independent 
of the scalar field operator in the matrix element. Provided that the operator does 
not introduce any significant exponential factor, the details of the operator will only 
influence the pre-factor, which will play only a small role in constraining the parameters 
of a given model. 

Physically, the importance of the virtual scalar field is that it provides a source of 
energy, making possible the transition to a critical bubble of true vacuum. In order to 



avoid further large suppression factors we should require that the virtual scalar field is 
able to efficiently supply energy to the region where the bubble is forming. This will 
be true if the region in which the energy of the virtual scalar field is concentrated is 
smaller than the size of the critical bubble. 



3 Width-Corrected Thin Wall Action 



In models of phenomenological interest the region of parameter space where false 
vacuum decay can be cosmologically significant is close to the thick wall limit p, 
1^, 1^. Therefore if we wish to apply the Landau- WKB method to cases of interest 
to phenomonology then we must extend the above calculation in order to take into 
account the effects of finite wall width. To do this we first need to derive a thin wall 
Lagrangian L{R, R) which includes the leading order corrections due to finite wall 
width. 

Following Voloshin ||^, we will consider a real scalar field with a potential given by 
Eq. (|T|). In the limit of small enough M and bubble radius much larger than the wall 
width, the thin-wall solution of the (p equation of motion may be written as 

1 



0(P) 



-V Tanh 



w 



-ip-k) 



:i8) 



where = — t^, w is the half width of the wall and k parameterizes the initial 
radius of the bubble. In order to calculate the width corrections to L{R, R) in the 
thin-wall limit (where R is now defined to be the radius at the point in the bubble 
wall at which = 0) we substitute the solution Eq. (18) into the scalar field action. 
This is given by 

(19) 



/ = -An J dt dr 
We write this a.s I = Ig + I^, where 



Ig = -47r j dt dr 



and 



-47r J dt dr (^-M^ 



(20) 



(21) 



On substituting the thin wall solution into these we obtain 



-271 J dt dp p{p^ + t^y/^Sech^ 
-AttM^v J dt dr r^Tanh 



1 



w 



(p-k) 



-ip-k) 
w 



(22) 



Evaluating these integrals then gives us a Lagrangian as a function of k and t, L{k, t). 
Since for the thin-wall solution Eq. ( ]T8| ) we know that the radius is given by 



R{t) = {k^ + t 



2.1/2 



(23) 



we see that k and t may be written in terms of R and R as 



k = R^l-R^ 



(24) 



and 



t = RR 



(25) 



Substituting these into L{k, t) then gives us the desired Lagrangian L{R, R) for which 



Eq. (23) corresponds to a minimum action solution. 

In order to obtain an analytic expression for the action from Eq. (^2]) we must 
fit the Sech'^(x) and Tanh(x) functions with approximate polynomial expansions. In 
practice we will use the expansion Sh4(a;) for Sech^(a;) and Th(x) for Tanh(x), where 
Sh4(a;) is defined by 



Sh4(a;) 



984 _ 1601 2 I _970_ 4 _ 1987 6 (I I ^ i n.\ 
1000 1000 ~'~ 1000 10000 ^\ \ I 



{\x\ > 1.5) , 



(26) 



and Th(a;) is defined by 



9937 



X 



Th{x) = 



10000 100000 

-1 {x < -2) 
1 (x>2) . 



28597 ^3 _|_ 62238 



1000000 



X° — 



571 ^7 
100000 



(|x| < 2) 



(27) 



In order to obtain the leading order width corrections, we expand the integrals in 
terms of ^. Such an expansion is possible for ^ small compared with 4, which is also 
necessary for the thin wall solution for Eq. (p!8D to be valid. 



Integrating Eq. (|22|), substituting Eqs. ^) for k and t and expanding in powers 
of w then gives the width corrected thin wall Lagrangian, 



L{R,R) = -jiyi-R^ 



+iR^+biRw^ (l - ^]+Oiw^) , (28) 



which is valid for \R\ ^ 4w and where a = 0.29 and b = 2.21. This reduces to the 
zero-width Lagrangian Eq. (|^) in the w ^ limit. 



4 Landau- WKB Calculation of the Catalysed De- 
cay Rate for Non-Zero Width 

In order to calculate the rate of catalysed vacuum decay we follow the same steps as 
for the zero-width calculation. The canonical momentum and Hamiltonian are now 
given by 



PR 



R' + cw^R^ 



- biRR^ 



(29) 



where c = 0.44, and 



?d3 



H + iR 



/i 



1 -i?2 



p2 

R^ + awHl - — + i?^) 



biR{l + ^)«;2 



(30) 



The parameters of the potential Eq. ([l|) are related to the parameters appearing in 
the Hamiltonian Eq. (^) by 



3 [32% 1 

Idtt 



, 3 \ 7 



3 flW^ 



= — 



47r/ {wfiy/"^ 



(31) 



We then solve the Hamiltonian equation in order to obtain i? as a function of R 
and E. In general we must solve 



E + beRR (l + Y J wA (1 - R^) = jl^ 



R^ + awHl -^ + R^) 



(32) 



This will give a 4th order polynomial m y = R^. To simplify this and so obtain an 
analytical expression, we will restrict attention to the case where the terms on 
either side can be neglected. In this case Eq. (|3^) reduces to 



/2l/' + /ll/+/0 = 



(33) 



where 

/2 = {beER + 2afl'^R'^)w^ 



/i = {E^ + biERw^ - ajl'^R'^ 
/o = -E^ + /i2^4 _ 2biERw^ + 2aj2^R^w^ . (34) 

Thus to this order one finds that 



(35) 



-fi + v/i - 4/2/0 

where 

/i -4/2/0 = E^ + {l3.26eE^R+11.7AE^jl^R^ -SMeEjl^R^ -2.32jl^R^)w^ 

+ {43.9Qi'^E^R^ + 3MeEfR^ - 1.26/i^/?^)w^ . (36) 

For sufficiently small width R^ can be expanded in terms of w"^ 

R? = Rl + w^Rj , (37) 

where i?Q is the zero width expression given in Eq. (||). In addition, so long as R^ is 
large compared with w'^Rf, we can expand Rq in powers of w"^, 

««fi„(l.lgi). (38) 

We will also assume that we can expand (1 — R^Y^'^ in powers of w"^. In this case we 
can give an analytic expression for the shift in the transition point x* relative to that 
obtained in the zero-width limit Xq = (— e/2)^^^. In terms of the rescaled variables 
e and d = we/fi, the canonical momentum is given by 

p{x, e) = po{x, e) (1 + (fg{x, e)^ — axRo(f , (39) 

where 



^(^,e) = i-,-A__ + 4^. (40) 



2 Rlil - Rl) ' x2 

The transition point x*(e) is then given by the solution of p{x, 0) = p{x, e). We define 
xi{e) by 

x*{e) =x*Q{e) + d^Xi{e) . (41) 



With the expansion p{x*, e) = po{x*, e) + (Pxip'q{x*, e), where 



Po(^,e) = ^-_=f , (42) 



(2v^e2 - x\ 

we obtain for the shift in the transition point 

, , 1.105 0.737 

The turning point x{E) is given by iJ(i? = 0) = E, so that Eq. (|30| ) imphes that 

x{Ef - x(E)2 + e + 6a;(E)d2 - ac/^ = , (44) 



which can be solved numerically with a and b as in Eq. 

To obtain the exponential suppression factor for catalysed decay, we evaluate the 
Landau integrals at the critical energy Cc so that in the scaled variables given by 
Eq. (0) 



eF(e„ d) = Re 



x*(ec) f^i^c) 

p(x, 0)(ix+ / p{x,ec)dx 



(45) 



10 Jx*{ec) 

with p{x,ec) in general given by Eq. (|39D. In the evaluation of the e=0 part of the 
Landau integral Eq. (^) there are two distinct regions. For > 1.15^^/^ 

we can 

expand and so p{x, 0) in powers of d"^ as in Eq. (ppj). For Ad ^ \x\ ^ 1.15d^^^ this is 
no longer possible and we find from Eq. (^) that R^ becomes 

R^^t^, (46) 

where is taken to equal +i in accordance with our definition of the phases in the 
complex plane. Thus the canonical momentum in this case is given by 

r,3 

p(x, 0) = ^ [{-ix"^ + 0.57xd) - 2.52i^/^x^/^d^/^\ . (47) 



We have integrated Eq. ( |i5| ) numerically, using the above expressions for p(x,e). 
The thin- wall approximation breaks down when |x| is smaller than d, and therefore 
we have used a small x cut-off at |x| = 4d in the rf-dependent part of p{x,e). This 
does not affect the final results appreciably. Some care must be taken in performing 
the complex integrals, in particular when the integration path is crossing branch cuts. 
Our result for the exponential suppression factor F{Ec, d), together with F{0, d) which 
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Figure 1: The scaled bounce action ^F{0, d) and the scaled action at the critical energy 
$,F{Ec, d) as functions of the scaled width d. 

corresponds to the bounce action, is shown in Figure 1. At (i = both curves agree 
with previous results |§. 

The quantity of particular interest to us here is the ratio of the exponential sup- 
pression factor at zero energy to that at the critical energy as a function of the wall 
width d. This measures the importance of catalysed decay, and we parameterize this 
by where 

F{Ec, d) = 0.16K(rf)F(0, d) . (48) 

In the limit of zero width /t(0) = 1 [||. 

We display K,{d) in Figure 2, and one can see that K{d) and the exponential sup- 
pression factor for catalysed decay increase as d increases, whilst the bounce action 
decreases. The increase in K,{d) with increasing d is in agreement with trend suggested 
by the numerical results for the thick wall bubble obtained by Kuznetsov and Tinyakov 
[0], which indicate a value K,{d) ^ 5 for the case of the thick-walled —X(j)'^ theory they 
considered. We have tested our results by comparing the d dependence of F(0, d) with 
that of the equivalent Euclidean bounce action calculated directly from the scalar field 
solution corresponding to a given value of d. We find that the Landau- WKB result 
at E = does indeed reproduce the bounce action for all d, up to the limits of our 
appr oximat ions . 
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Figure 2: The parameter k as a function of d. 

In obtaining these results we have made a number of approximations, which im- 
pose hmits on the range of wall widths we are able to discuss. We find that the 
tightest constraint comes from the requirement that the expansion of in powers of 
w, Eq. (|38|), is valid at the transition point, which requires, for the e = part of the 
Landau integral, that 2.32d^ < |a;*(ec)|'' for |x*(ec)| ~ |a;o(<2c)| = 0.43. This requires 
that d is small compared with 0.05. Although this is a tight constraint on the range 
of widths for which our calculation is valid, physically we do not expect any change 
in the d dependence of n{d) at (i ~ 0.05. We should be able to extrapolate the d 
dependence of n{d) we obtain at small values of d, for which our calculation is valid, 
up to larger values of d so long as the wall width is small compared with the radius 
of the bubble, which provides the only other length scale in the problem. Fitting hi{d) 
with a polynomial up to d^, using values of d up to 0.01, we obtain 

K{d) ^ 0.999 + 0.125rf + llO.Srf^ + (-4652^^ + 135656rf^). (49) 

The d^ and c?^ terms only become important once d ^ 0.02, indicating that they are 
artifacts of the breakdown of our approximations. Therefore we see that the dominant 
behaviour of K,{d) at small d is d"^. Thus we can extrapolate K,{d) in the form 

K{d) ^ 1 + llOd^. (50) 



5 Cosmic Ray Catalysed vs. Spontaneous False 
Vacuum Decay 

In this section we give the condition for cosmic ray catalysed vacuum decay to be more 
important than spontaneous false vacuum decay. Cosmic rays can catalyse electroweak 
vacuum decay by scattering off background protons and producing a virtual Higgs field 
of sufficient energy (which will be spherically symmetric in the centre of mass frame). 
As discussed earlier, the importance of the virtual scalar field is that it provides a 
source of energy for the growth of the true vacuum bubbles. For the virtual scalar 
field to be able to provide all of its energy to the critical bubble without introducing 
additional large suppression factors, it should be concentrated in a region at least 
as small as the radius of the critical bubble. In fact, for the thin-wall bubbles we 
are discussing here, the Compton wavelength of the virtual scalar field at the critical 
energy turns out to be generally less than or of the order of the width of the bubble 
wall, ^ w. This follows because the condition ^ w is equivalent to 



where Vf, = Xv^/A is the height of the potential barrier. Thus with Cc ~ 0.15 this will 
be satisfied if e ^ 3.7Vb/X. Since e is the splitting of the energy density of the vacuum 
states, which is less than or of the order of Vf, for the case of thin-wall bubbles, we see 
that this will generally be satisfied in the thin-wall limit we are considering. 

The decay rate is largest for cosmic-ray collisions of centre of mass energy ap- 
proximately equal to the critical bubble energy. This is because, as discussed in the 
previous sections, the scalar field catalysed decay probability is highest for a scalar 
field of energy equal to the critical energy E^. Since the fiux of cosmic rays rapidly 
decreases with increasing energy |Ty] whilst the probability of vacuum decay does not 
increase with energy once the scalar field energy is greater than the critical energy 0], 
it follows that the contribution to the vacuum decay rate from higher energy cosmic 
ray collisons will be negligible. 

The cross-section for a proton-proton collision to produce a Higgs boson of energy 




(51) 



Eh via gluon fusion is given by [11 



'^'^ {E,^, Eh) = , ^ ^J'^l^ x,X2G{x,)G{x2) , (52) 



where G(x) is the gluon distribution function, xi = (-^)e^, X2 = (-^)e~^ and 



e^y = i^Eh + y El~ ml) {Eh — y E^ — ml). Thus if we are interested in producing a 
Higgs boson in an energy range 6 Eh around Eh for a given Ecm, the cross-section is 
given by 



aiE,^,Eh) ^ -^SEh , (53) 
dEh 



which gives, for E^ ^ ml, 



2 

where Xi = |^ and X2 = tt^-w — • We will consider SEh ~ Eh in the following. The 
centre of mass energy is Ecm = y^2mp£'i, where Ei is the cosmic ray energy and mp 
the proton mass. Since it is best to produce a Higgs of energy Eh with the minimum 
cosmic ray energy Ei, we shall consider the cross-section for the case E^m ~ Eh, 
corresponding to cosmic rays of energy Ei ^ • The total number of cosmic ray 
induced vacuum decay events during the history of the Universe is then given by 

N{Eh) ^ f{E{)a{E{)P{Eh)n,ptJE, , (55) 

where f{E) is the flux of cosmic ray protons of energy E, P{Eh) is the probability 
of catalysed vacuum decay in a virtual scalar field of energy Eh, Up ^ 2 x lO''^ is the 
number of protons in the Universe and ^ 5 x 10^^ GeV^^ is the age of the Universe 
[0. In this we have considered N{Eh) to come mostly from a range 6E1 ^ Ei around 
the minimum cosmic ray energy needed to produce a Higgs of energy Eh, this is 
consistent with 6 Eh ~ Eh- The cosmic ray flux is given by (for Ei > 100 GeV) |1C 



E 



-7 



where A = 2.2 x 10"^^ GeVl We define x{Eh) by N{Eh) = x{Eh)P{Eh). From 
Eq. (|55|) we find 

x{Eh) ~ e^^^^-''^^'' , (57) 



where 

Thus, with Tj defined by P{Eh) = e~^, we find that catalysed false vacuum decay will 
typically occur for r] ^ 130. 

Spontaneous false vacuum decay occurs at a rate per unit volume per unit time 
given by H 

^spont ~ ml exp (-5'4) , (59) 

where (= F{0,d)) is the bounce action. Thus the probability of spontaneous 
vacuum decay during the history of the Universe is given by 

Tspontturl ^ lO^^^e-^^ , (60) 

where ^ tu is the radius of the Universe and we have used rrih = 10^ GeV. Therefore 
spontaneous vacuum decay will occur if 5*4 ^ 403. In terms of n{d), since F{0,d) 
is equal to the bounce action 5*4 and since F{Ec, d) ^ 130 implies that cosmic ray 
catalysed false vacuum decay occurs, we see that catalysed false vacuum decay will be 
more important than spontaneous decay if 

K{d) < 2.0 . (61) 

Using the extrapolation of K{d), Eq. (0), we find that K{d) ^ 2 occurs at d ^ 0.1, 
well outside the thick-wall limit where our results would become unreliable. (The 
initial radius of a bubble at the critical energy is Xc = 0.67 whilst the total thickness 
of the wall is 2d. Therefore the bubble wall width in this case is about 30% of its 
initial radius). Thus we can conclude that cosmic-ray catalysed false vacuum decay 
can only be important if the bubbles are not thick-walled. In terms of the parameters 
of the potential Eq. (|l]), the condition Eq. (|6TD implies that ^ O.lAf^. 

Could false vacuum decay to thin-walled bubbles ever be cosmologically significant? 
The answer is almost certainly not, at least not for the case of a perturbative model. 
To see this, consider spontaneous vacuum decay. For a real scalar field with a potential 
given by 

Vi<t>) = - 1<I>' + (62) 



to which most phenomenological models can be roughly fitted, the bounce action is 
given by ^ 

S, = l^S,{X), (63) 

where 5'4(A) is the rescaled bounce action as a function of A and A ~ |m^A/5^. Thin- 
wall bubbles correspond to A ^ 0.85, for which d ^ 0.2. In this case it may be shown 
that ^4(A) > 2600. Therefore 

^4 - ^ > (64) 

where we have used 5^ ~ |m^A for the thin-wall limit. Since cosmologically significant 
spontaneous false vacuum decay requires that 5*4 ^ 400, we see that we would require 
A ^ 3.3 in the thin-wall case. Since for the theory described by Eq. (^) the 1-loop 
correction to the scalar self-interaction becomes larger than the tree-level contribution 
for A ^ 2, we see that thin-wall false vacuum decay could only be cosmologically 
significant if the scalar field were non-perturbative. Although this result is true for 
spontaneous false vacuum decay, the possible enhancement of the cosmic ray catalysed 
false vacuum decay rate over the spontaneous decay rate in the thin wall limit will not 
alter this conclusion, since the increase in 5*4 (A) is far greater than the reduction in 
K{d) as we go to thinner bubbles. 



6 Conclusions 

In this paper we have attempted to address the question of whether false vacuum 
decay catalysed by cosmic ray collisions could ever be an important constraint on 
models of particle physics. We have used the Landau- WKB method, including the 
effects of finite bubble wall width, to calculate the rate of cosmic ray catalysed false 
vacuum decay relative to that for spontaneous decay. We find that cosmic ray catalysed 
decay is the dominant process for sufficiently thin-walled bubbles. However, for the 
case of perturbative models with cosmologically significant false vacuum decay rates, 
which almost certainly decay via thick-walled bubbles, our results strongly suggest 
that spontaneous false vacuum decay via tunnelling will always be the dominant decay 
mode. 



We should also note that our results indicate that vacuum decay is very unlikely to 
be catalysed by particle collisions in accelerators; the total luminosities involved are 
simply far too low. In the very early universe the situation might however be differ- 
ent. It is conceivable that, for a certain range of the potential parameters, catalysed 
decay could be more important than thermal fluctuation over the barrier or the zero 
temperature tunnelling. Such a possibility could constrain particle physics models in 
an interesting way and merits further study. 

Our results are based on extrapolating the cosmic ray catalysed decay rate from the 
small values of wall thickness for which our calculation is valid to larger thicknesses 
(although still well outside the thick- wall limit). This extrapolation indicates that 
catalysed false vacuum decay becomes less important than spontaneous decay once 
the bubbles produced by cosmic ray collisions have a thickness greater than about 
30% of their initial radius. It is possible that the behaviour of the catalysed decay rate 
as a function of wall thickness could begin to deviate from the thin-wall behaviour at 
such thicknesses, although we do not expect our conclusions to change significantly. 
However, it would be valuable to have a more accurate calculation for wall thickness 
which are between the thin and thick wall regimes. To go beyond the thin wall 
approximation using the Landau- WKB method would require an action describing 
the dynamics of thick(er)-walled bubbles as a function of a single variable R. This 
could be done, for example, by considering numerical solutions of the (p equations 
of motion and then deducing an Ansatz for (p which is a reasonable approximation 
to the observed behaviour. Substituting this into the scalar field action would, in 
principle, allow for the derivation of a thick wall action as a function of R, to which 
the Landau- WKB method could be applied. 

The analytical Landau- WKB approach we have considered here is complementary 
to the numerical approach as followed by Kuznetsov and Tinyakov [0]. The trend of our 
results with increasing wall width agrees with that suggested by their thick wall bubble 
results. We believe that the Landau- WKB approach can give us additional insights 
into the physics of catalysed false-vacuum decay and provides a powerful approach to 
the problem. 
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Appendix. Landau- WKB Method in Quantum Me- 
chanics 



In his original discussion, Landau was concerned with calculating the matrix ele- 
ment of an operator f{x) between two states, ipi and tp2, corresponding to particles 
of energies Ei and E2 {E2 > Ei) reflecting from a potential barrier U{x). The pres- 
ence of the potential barrier makes it possible for the matrix element between the two 
energy eigenstates to be non-zero. The standard WKB wavefunctions for this system 
are given by (with h = 1) 



^1 

and ip2 = i't + "^2" where 

i^t = 



exp 



Pi dx 



cos 



/ Pidx-j 

Jai 4 



X < ai 



X > tti 



(65) 



-iCo 



y\p2\ 
C2 



exp 



exp 



P2 dx 



a2 



X < a2 



/■^ , in' 



X > a2 



(66) 



and ip"! = i^'t)*- III this the full 1^2 wavefunction has been set to zero at x < 02 < oi, 
where it does not contribute significantly to the matrix element (both ipi and 1IJ2 being 
exponentially suppressed in this region). The matrix element of an operator / between 
ipi and 4'2 is then given by /12 = /12 + /12, where 



12 



'ipifi't dx 



(67) 



and /12 = (/i2)*- The WKB wavefunctions are not valid close to the turning points ai 
and a2. However, by continuing x to the upper half of the complex plane these points 
may be avoided. The wavefunctions in the upper half-plane may then be written as 



^1 



2[ipiYl^ 

-1C2 

2[ip2V/' 



exp 
exp 



X 

I I pidx 

ai 

X 

—% I p2dx 

0.2 



(6^ 



The matrix element is then given by 

-lCiC2 



fl2 — 



f{x)dx 



X exp 



i I pidx —if P2dx 



(69) 



By the residue theorem, the matrix element will be determined by the points in the 
complex plane ("transition points") where the true wavefunctions are singular, which 
corresponds to the points at which the exponential factor is stationary as a function of 
X, which we denote by x*. This follows because at these points the canonical momenta 
for two different energies arc equal, which implies that the momenta are independenent 
of the energies of the states. Therefore the potential and so the wavefunctions at these 
points must be singular. The integral, up to a relatively unimportant prefactor, will 
then simply equal the exponential factor evaluated at x*. Thus 



12 



exp 



r ^2m(C/ - El) dx - r ^j2m{U - E2) 



dx 



(70) 



where we have used Pi = —iy2m{U — Ei) and where the square roots are taken to be 
positive on the real axis for x < a2. Since the complex part of the exponent contributes 
only an insignificant phase factor, we can write /12 = /12 + /12 ^ 



fi2 ^ exp 



Re 



2m{U - El) dx- ^j2m{U - E2) dx 



(71) 



This expression is true so long as the states in the matrix element are "strongly dif- 
ferent" , meaning that the exponential suppression factor is large. 
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